Spin and orbital effects in a 2D electron gas in a random magnetic field. 
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Using the method of superbosonization we consider a model of a random magnetic field (RMF) 
acting on both orbital motion and spin of electrons in two dimensions. The method is based on 
exact integration over one particle degrees of freedom and reduction of the problem to a functional 
integral over supermatrices Q(r, r'). We consider a general case when both the direction of the RMF 
and the g-factor of the Zeeman sphtting are arbitrary. Integrating out fast variations of Q we come 
to a standard collisional unitary non-linear a-model. The collision term consists of orbital, spin and 
effective spin-orbital parts. For a particular problem of a fixed direction of RMF, we show that 
additional soft excitations identified with spin modes should appear. Considering 5-correlated weak 
RMF and putting g = 2 we find the transport time rtr- This time is 2 times smaller than that for 
spinless particles. 
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I. INTRODUCTION 

Models of a random magnetic field (RMF) acting on 
electrons in two dimensions are intensively studied in 
mesoscopic physics. There are direct experiments on 
high-mobility heterostructures subjected to a magnetic 
field of randomly pinned flux vortices in a type-II su- 
perconducting gate^, type-I superconducting grains^ or a 
demagnetized ferromagnet"^. For theoreticians, the RMF 
models are important as an example of systems with an 
interaction reduced to a gauge field. These models arise 
in theory of composite fermions in the fractional quantum 
Hall effect near half-filling'* as well as in a description of 
doped Mott insulators^. 

From the theoretical point of view, one of the most in- 
teresting tasks in the study of any model with a disorder 
is to determine the large scale behavior of electrons and 
find universal properties of the model. In doing so, one 
may come either to a metal or insulating behavior. In 
the latter case electron wave finctions should be local- 
ized. In the context of the RMF models the localization 
have been discussed in many numerical works, which re- 
sulted in three different conclusions: a) localization of 
all the states^"^; b) existence of a band of delocalized 
states^ ~^^; c) localization of all the states except those 
in the band center*^ ^*^. 

Analytically, the RMF models were studied by the 
standard diagrammatic technique*^ as well as using dif- 
ferent non-linear cr-models^^~^*. The final conclusion 
drawn by using the methods was that the RMF models 
belonged to the usual unitary class of universality with 
localization in two dimensions (provided the correlations 
of the RMF {BqB-q) at small g ^ increase slower than 
1/q^ , Ref.20). 

Ususally, when deriving the proper a-model for the 
RMF problems one uses the standard scheme^^ based on 
the saddle-point approximation. Calculations presented 



in Refs.^®"^*^ are performed in this way. From the point 
of view of the conventional perturbation theory this ap- 
proximation corresponds to the self-consistent Born ap- 
proximation (SCBA), which is not good for a long range 
disorder. The same approximation has been used in the 
diagrammatic approach of Ref.^*. 

In addition to the difficulty of a description of the long 
range disorder, the standard scheme is not convenient 
for a generalization of the RMF model for the case when 
spin degrees of freedom are important. This is because 
the effect of the magnetic field on the orbital electron mo- 
tion is accounted for by adding a vector potential in the 
Hamiltonian, whereas the interaction with the electron 
spin is described by the magnetic field itself. As the cor- 
relations of the vector potentials and the magnetic fields 
are different, it is not easy to consider both the effects on 
equal footing. 

Partly this problem has been resolved in a recent 
paper^^ for free electrons in 2D with the g-factor g = 2 
and a magnetic field perpendicular to the plane. Un- 
fortunately, a mathematical trick of replacing the initial 
Hamiltonian by a Dirac Hamiltonian used in that paper, 
which was the basis of the suggested calculation scheme, 
can be applied only for this particular system. 

Another analytical method suggested in Ref.^^ enables 
one to avoid using the saddle-point approximation and 
to obtain a ballistic a- model applicable at all distances 
down to the Fermi wavelength Xp. This method is based 
on using quasiclassical equations of motions which con- 
tain not the vector potential but the magnetic field only. 
As the vector potential itself does not enter the ballsitic 
cr-model one can include rather easily the Zeeman term 
without extra assumptions about the value of the ^-factor 
and the direction of the magnetic field. However, this is 
still not the most general method because it is essentially 
based on the quasiclassical approximation and therefore 
short range correlations of the magnetic field cannot be 



considered. 

In the present paper we use a new method of super- 
bosonization suggested recently^'*. The method is based 
on the exact integration over the one particle motion and 
reformulation of the initial fermionic problem in terms of 
a functional integral over supermatrices Q{r,r'). As it 
has been shown in Ref.^**, in the quasiclassical regime 
(smooth disorder and lengths larger than Xp) the ma- 
trix Q{r,r') corresponds to the matrix (5n(r) of Ref.^'"'. 
The method of the superbosonization has several advan- 
tages. First, it uses neither saddle-point nor quasiclas- 
sical approximations and is exact. Both short and long 
range disorder can be considered on its basis. Second, 
integration over one particle motion is carried out before 
disorder averaging. The latter enables one to consider 
interaction effects related to the random scattering more 
carefully. Finally, the method results in a simple expres- 
sion for the energy in terms of the matrices Q(r, r') and 
makes the disorder averaging a rather simple procedure. 

The new theory is invariant with respect to rotations 
of the matrix Q{r,r') in the superspace provided they 
commute with the Hamiltonian H . This makes possible 
separation of the massive modes from the soft ones and 
we can integrate them out. We show that this proce- 
dure is justified only in the regime when the scattering 
effects do not result in a strong coupling between elec- 
trons. For the short range disorder this is so over dis- 
tances exceeding the correlation length of the disorder 
(or the Fermi wavelength Xp) whereas for the smooth 
disorder the coupling becomes sufficiently weak beyond 
the Lyapunov length. The latter case was discussed in 
many details in Refs.^^'^^. 

The low energy theory found is described by a ballis- 
tic nonlinear a- model with a collision term. This term 
consists of three parts that can be related to the orbital, 
spin and some effective spin-orbital scattering. Similarly 
to the model with magnetic impurities (see Ref.^^), the 
second scattering results generally in the relaxation of 
all spin modes. At the same time, we show that in the 
model with a fixed direction of RMF an additional soft 
mode corresponding to fluctuations of the spin along the 
field should appear. Finally, using standard method of 
integrating out the angle modes we come to the conven- 
tional unitary diffusive a- model and find the transport 
time Ttr- 

The paper is organized as follows. In Sec. II we discuss 
the superbosonization procedure and give an alternative 
derivation of supcrbosonized theory. At the end of the 
section we discuss symmetry properties of the obtained 
model and derivartion of the non-linear a- model. 

In Sec. Ill we introduce the main definitions of RMF 
model involved, derive the a- model valid in the coUi- 
sional regime and discuss conditions of its applicability. 

In Sec. IV the same problem is considered in the dif- 
fusive limit. We calculate with the help of the obtained 
diffusive cr-modcl the transport time and the spin suscep- 
tibility. 



II. SUPERBOSONIZATION AND DERIVATION 
OF THE CT-MODEL. 

Below we consider a two-dimensional {2D) electron gas 
placed in a static inhomogeneous magnetic field. The 
field is assumed to act both on the orbital motion and 
electron spin. Our consideration is based on the new 
method of superbosonization proposed in recent paper^**. 
This method uses exact integration over electron degrees 
of freedoms and reformulation of the problem in terms of 
integrals over supermatrices with a rotational symmetry 
(nonlinear cr-model). 

Before starting the study of RMF problem we would 
like to give alternative derivation of the supcrbosonized 
model. Although, the scheme of the derivation presented 
in Ref.^** is straightforward and exact, the final repre- 
sentation of the Green functions in terms of a functional 
integral over the supermatrices can be obtained even in 
a more simple way. Following the standard way (see the 
book^^) we introduce first a generating functional Z[a]: 



Z[a] = / exp(-La[^])i?V 



(2.1) 



where the Lagrangian La [V'] has the form 



LM = -I / V'(r) ( Fr - e + I + ^^A ) ^(r)dr + 



i I V'(r)a(r,r )i/;(r )dvdY 



(2.2) 



Hr is Hamiltonian of the initial model (we write it below) 
and ■)/'('") is a supervector. The conjugated supervector 
VJ(r) is related to ^(r) according to the following defini- 
tion: 



i,{v) = (CV')^(r) 



(2.3) 



The structure of the supervectors '0(r), "0(1") as well as the 
form of the matrix C depends on the problem involved 
(for details see Ref.^^). For the problem of electrons with 
spin considered below ^(r) should be a 16-component su- 
pervector. 

The Lagrangian LaiV'li Eq.(2.2), differs from the con- 
ventional one by the presence of the source a(r, r'). In 
general, the source term is for the problem considered an 
arbitrary 16 x 16 supermatrix depending on two coordi- 
nates r, r'. Differentiating in elements of this matrix one 
can obtain various correlation functions. For example, 
the level-level correlation function R{uj) can be written 
as 



R{u 



1 



1 



lim Re- 



92 



2(7ri^F)2 ai, 0:2=0 daida 



-Z[a] (2.4) 



provided the source a(r, r') is taken in the following form 
a(r, r') = a(r)(5(r — r') 



a(r) 



di 
—6(2 



ai,2-^(l-fc), (2.5) 



fc = ±1 in fermionic and bosonic blocks respectively. As 
supermatrix ^Q(r)V'^(r') is self-conjugated one may con- 
sider only self-conjugated sources: 



a(r,r') = Ca^(r',r)C^ = a(r,r') 



(2.6) 



Sources with constraint Eq.(2.6) have minimal number of 
elements required when finding an arbitrary correlation 
function. 

The integral defined in Eqs.(2.1), (2.2) is gaussian and 
can be readily calculated: 



Z[a\ — exp ( — Strln 



H - e-\ 1 A - a 

2 2 



(2.7) 



The logarithmic derivative of the partition function Z[a\ 
Eq.(2.7), in a is a matrix Green function G(r, r'): 






satisfying the following equation 

LU LU + iS 



(2.8) 



Hr-e 



A- a] G(r,r') = i<5(r-r') 



(2.9) 



Calculating the Green function G(r, r') from Eq. (2.9) 
and using Eq.(2.8) one can find the partition function 
Z[a], Eqs.(2.1),(2.7). 

What we want to show now is that the Green function, 
Eq. (2.9), can be represented exactly as an integral over 
supermatrices Q (r, r') in the following form 

G{r,r') ^ Z-'[a] / 0(r, r') exp(-F„[Q])i?Q (2.10) 



(2.11) 



and the functional -FaiQ] has the form 



where Z[a\ is a new partition function 
Z[a]^ f c^p{-Fa[Q])DQ 



S{r-r')Q{r,r')drdr'+ (2.12) 



-Str In Q - -Str / a(r, r')Q(r', r)drdr' 
z z J 

Supermatrix Q(r, r') in the integral Eq.(2.10) is not arbi- 
trary and should have a certain structure. First, it must 
be self-conjugated: 



Q(r,r') = g(r,r') 



(2.13) 



to provide the same for the result of the integration in 
Eq.(2.10). 

Second, the structure of the supermatrix Q{r,r') 
should be defined such that the energy Fa[Q], Eq.(2.12), 
would have a minimum. One can show that this takes 
place if Str((5^) is non-negative. The latter requirement 
is fulfilled for the supermatrices with the following con- 
straint: 



Q^{r,r')^KQ+(r\v)K, K 



1 
k 



(2.14) 



where Q± = (1/2)A[A, Q], k is the same as in Eq.(2.5). 
Relation (2.14) will be used below when formulating the 
structure of a cr model. 

Eqs. (2.10-2.12) have been suggested in Rcf.^^ and we 
want to demonstrate now a simpler way of their deriva- 
tion. In order to prove Eq.(2.10) we write the following 
identity 



-2iZ-'^\a] 



(5exp(-iStrlnQ) 
<5Q(r",r) 



Q(r",r')dr" 



exp 1 -2^*^^ 



LO LO + iS ^ 

H -e^ \ A - a 

2 2 



Q]DQ^ 



= iS{r — r'). 



(2.15) 



and integrate over Q by parts. The derivative S/SQ 
should act now on both Q and the exponential. At this 
point, the supersymmetry plays a crucial role. Differenti- 
ating first the supermatrix Q we obtain the supermatrix 
product (S/SQ) Q. As the number of the anticommuting 
variables in the sum over the matrix elements is equal to 
the number of the boson ones and the derivatives have 
the oposite signs, this matrix product vanishes. Differ- 
entiating the exponential only we come to the following 
equation: 

Z~y] I dv" fH-e+^ + ^^A - a) (r,r")x 



Q(r", r') exp(-F,[Q])7^0 = iS{r - r' 



(2.16) 



Eq. (2.16) proves immediately that the integral Eq.(2.10) 
does satisfy Eq.(2.9) and we have really the alternative 
representation of the Green function in terms of an inte- 
gral over the supermatrices Q. 

Integrating Eq.(2.10) over the source a(r, r') we con- 
clude that partition functions Z[a], Eqs. (2.1, 2.7) and 
Z[a], Eq.(2.11), are proportional to each other with a 
coefficient that is independent of the source. Putting 
a(r,r') = in the definitions of the both functions and 
finding that Z[a = 0] = Z[a = 0] = 1 due to supersym- 
metry we come to the equality: 



Z[a] = Z[a] 



(2.17) 



We emphasize that the relation (2.17) is exact and does 
not depend on the form of the Hamihonian Hj-, Eq.(2.2), 
and the source a(r, r') provided both the partition func- 
tions are defined in Eqs.(2.1), (2.11) in terms of the su- 
perintegrals. 

In the low energy limit, the field theory specified by 
Eq.(2.12) can be reduced to an effective nonlinear a- 
modeP'*. A possibility of this reduction is related to the 
invariancc of the free energy functional Fa with respect 
to the rotations of the supermatrix Q 



Q^VQV 



(2.18) 



for the case when the unitary matrix VV = 1 commutes 
with the Hamiltonian iJ, Eq.(2.2). This means that the 
low energy limit of the model is determined by rotations 
V with small values of commutator [iJ, F]. For differ- 
ent problems the smallness of the commutator can be 
provided by imposing on the rotations V different condi- 
tions. For RMF model we discuss this point at the end 
of the section III. 

In order to reduce the general formula, Eq.(2.12), to 
a a- model containing only the rotations V we represent 
the supermatrix Q in the form 



VqV, VV = 1 



(2.19) 



where g is a diagonal matrix [q^ A] = 0. One can impose 
the condition Vh = KV to fix the gauge freedom re- 
lated to the invariance of the definition, Eq.(2.19), with 
respect to the replacement V -^ Vv, where [f,A] = 0. 
The spectrum of fluctuations of the matrix q has a gap 
and they can be considered using the saddle-point ap- 
proximation. Substitution of Eq.(2.19) into the free en- 
ergy functional, Eq.(2.12), and integration over q carried 
out in the quadratic in fluctuations of q approximation 
results in the following expression for the effective free 
energy functional 



F[V] = -Str ( V 



'h + ^^a,v 



gl-ht4V[H,V]gY, 



(2.20) 



where g is the Green function, Eq. (2.9), in the absence of 
the source a — 0. This expression has been first obtained 
in Ref.^^ and is just the first two terms of the expansion of 
the free energy functional in powers of the commutator 
[iJ, T^]. As it has been shown in Ref.^**, this approxi- 
mation is sufficient if the scattering is rather weak. We 
use Eq.(2.20) as the starting point for the study of RMF 
problem and consider the limit of a weak field. 



III. RMF: FORMULATION OF THE PROBLEM. 
REDUCTION TO COLLISIONAL a-MODEL. 



in Eqs.(2.1), (2.2) 16-component supervectors ip{r). The 
most convinient choice for their structure is as follows 
(see also^^): 



^ 






r 



1 



X'" 



1 

71 



iCFyS'''' 



(3.1) 



where m — 1,2 divides the supervector space into ad- 
vanced (A) and retarded (R) subspaces; x™j 5"" are anti- 
and commuting 2-component vectors in spin space, Uy is 
the second Pauli matrix. The matrix C that determines 
the charge conjugation is written now as 



C = A. 



ci 

C2 



Cl 



C2 = 







(3.2) 



where the matrix A is the third Pauli matrix in the 
advanced-retarded space. Below we use also the matrix 
T3 = ±1 in the time- reversal space and 11 = t^ ® a, 
c — (ca:, CTy, cr^). Using these definitions we write the 
Hamiltonian H^, Eq.(2.2), as follows 



H,= 



1 
2m 



-iVr - -T3A(r) 
c 



- e^ - |AiBB(r)S, (3.3) 



fiB = e/{2mc), g is a factor that determines the Zeeman 
splitting. In Rcf.^'' the authors put g = 2 and assumed 
that the magnetic field B(r) was perpendicular to the 
plane of electron gas. This was the only case when the 
mathematical trick used in that paper and based on re- 
placing the initial Hamiltonian (3.3) by the Dirac Hamil- 
tonian could be realized. In the present paper, using 
the new scheme of the calculations we consider a more 
general case with an arbitrary value of ^-factor and the 
direction of the magnetic field. 

As mentioned previously, we consider the limit of a 
weak magnetic field. This assumption implies that the 
corresponding radius Rh of the cyclotron motion Rh = 
vp{eB /rac)^^ , vp is Fermi velocity, is larger than other 
characteristic lengths. Then, considering electrons in a 
region with sizes smaller than R^ one can choose the vec- 
tor potential A(r) such that it would be small everywhere 
in this region. The smallness of the vector potential al- 
lows one to represent the Hamiltonian Eq.(3.3) as the 
sum of the main part Hq and a small perturbation 6H 



In order to take into account the spin of the electrons 
one should double the number of variables and consider 



Hr — H( 



Or 



6Hr 



Hor — —- 



2m 



EF, 



SHr 



^£^^^(^)^r+2^^'W-f^-^^W (^-4) 



We choose for the vector potential the following gauge: 
divrA(r,z = 0) = d^A^{r, z = 0) + dyAy{r, z = 0) = 0, 



gop is the Fourrier transformed Green function Eq.(2.9) 
in the absence of the field, A = and the source a = 0. 
We find the commutator [7Yo,V^] and write it in the 
Wigner representation as 



t[no,v] = — yp(r) + — - — 



[A,Vp{r)] (3.10) 



A,{r,z = 0) = 0, 

r is coordinate on the plane of electron gas. This gauge 
corresponds to the well-known London gauge in the su- 
perconductivity theory. 

The next step in the calculation is to average the par- 
tition function Z[a\ = J exp{-F[V])DV with F[V] from 
Eq.(2.20) over the random magnetic field. In the present 
paper we take a gaussian distribution of the field with 
the pair correlation 



{Bdr)B,(r'))=2(—Yw,,{r-r') 



(3.5) 



The free energy -F[y], Eq.(2.20), is written in a somewhat 
implicit way. In order to reduce it to a more explicit way 
we use the Wigner representation. 
For any matrix O it is defined as: 



0(r + p/2, r' - p/2) = J Op{r)e'PPdp 



(3.6) 



The product of two operators Oi, O2 can be written in 
the Wigner representation in the following form: 

Oip(r)*Q2p(r) = 



Oip(r)exp 



^r'^p - ^p^i 



02p(r) (3.7) 



The new operation * introduced in Eq.(3.7) has all prop- 
erties of the usual matrix product. In particular, it is 
associative. 

Using the Wigner representation we expand Eq.(2.20) 
in 5H and regard only terms up to the second order 



F[V]^F„[V]+F^[V]+F2[V], 



(3.8) 



where r.h.s. is the sum of terms of the zeroth, first and 
second order in the field, respectively. Below we present 
an explicit calculation of the zeroth term only. Calcula- 
tion of the other terms can be carried out in the same 
way. 

The zeroth order term Fg \V] takes the form 

MV] = '-StT{V[Ho,V]go) - ^Str{V[Ho,V]go)\ 



p2 , UJ + lS ^ 



(3.9) 



Eq.(3.10) contains a small frequency lo and the gradi- 
ent VrVp(r). Being interested in variations at distances 
much exceeding the wave lentgh we neglect the second 
term in Eq.(3.9). In the same approximation one may 
replace the "star product" * Eq.(3.7) by the usual one 
and write the first term in Eq.(3.9) as follows 



Fo[V] = i 



Vp(r) 



pVr i{uj + iS) 



A 



m 



Vpir)gopdrdp, 
(3.11) 



The Green function gop has a sharp peak at the Fermi 
surface, whereas the function Vp is smooth. Therefore, 
one may replace in Eq.(3.11) the momentum p by its 
value at the surface: p -^ Pfh, consider separately the 
integration over the absolute value |p|, which is equiva- 
lent to integration over ^ = p^/2to — ep, and unit vector 
n = p/p. Making the replacement dp -^ lyd^dn and 
carrying out the integration over S, only in the Green 
function gop we have 



9opd^ 



^-e+-{l + A)+iSA dS,=TrA 

(3.12) 



Using Eq.(3.12) we come to the following result for Fq, 
Eq.(3.11) 



Fo[V] = ^Str 



AF„(r) vpnV 



i{uj + i5) 



A ) Vn{v)dvdTa., 
(3.13) 



This expression is just the usual kinetic term of the bal- 
listic cr-model with the matrix (5„ ~ V-aAV^i- In order to 
take into consideration the scattering one should calcu- 
late the rest terms Fi[V], F2[V] of Eq.(3.8). 

In the main in small commutator [Hq^V]^ Eq.(3.10), 
approximation the term i^i[y], Eq.(3.8), may be written 
in the following form: 



Fi[V] = -SiT{V[5H,V]g^) 



(3.14) 



Using the Wigner representation and neglecting high gra- 
dients of the matrices Vp (r) , Vp (r) we find with the help 
ofEq. (3.4) 

FAV] = i^Str j e^^^ [-r3pA(q) - (5/4)SB(q)] x 



^p-§ (r)3op Vp+§ (r)c?rrfpdq, 



(3.15) 



In Eq.(3.15), we disregard the term (e^/2TOc^)A^(r) in 
5Hr, Eq.(3.4), because it results in small in the magnetic 
field corrections. Moreover, considering large distances 
and therefore assuming that the matrices V-p{v) vary in 
space slower than the vector potential A(r) and magnetic 
field B(r) we conclude that the functional -Fi[T^] should 
be self-averaging and vanish. Thus, only term F2 \V] in 
the expansion Eq.(3.8) gives a contribution related to the 
scattering in the RMF. In the main order it reads: 

F2[V] = -^StT{V[SH,V]go5Hgo)-^StT{V[SH,V]goy « 



i^y^^'J [-^3PA(q) + (g/4)SB(q) 



^p-f (r).9op-f Vp+t« hr3pA(q') + (.9/4)EB(q')] 



V^_^{T)g^^_^V^^^ir)e^'^^+^'^''drdpdqdq' (3.16) 

Eq.(3.16) can be easily averaged over the magnetic field. 
Noticing again that gop is a sharp function of the mo- 
mentum p at the Fermi surface and using Eq.(3.5) we 
come to the a model with the following effective energy: 

-^ [^JnJ — ^ kin[^n\ i ^ orb[^n\ i ^ sp[^n\ ' ^ sp—orb[^n\7 

where Fkin [Qn] determines ballistic (or free) propagation 
and coincides with i^o[Qn] Eq.(3.13) whereas the other 
terms are related to three types of scattering which can 
be called as orbital, spin and spin-orbital, respectively: 



ForbiQn] = z — ^*'' / '^ ^^Wzz(ni - n2) 



8 7 1 — 111 n2 

X Qni (r)Qn2 {r)drdnidn2, 



(3.17) 



^ sp [^/r 



{Triyf 



Str 



ujy(ni -n2) X 



Qm {r)<TiQn2 {r)(Tjdrdnidn2, 



(3.18) 



-*■ sp—orb [ 



(TTiy) 



-Str 



g\ [ni X n2]. 



2 7 V 4/ 1 - nin2 

X crjQni(r)(5„2(r)drdnidn2, 



(3.19) 



where Wij (ni — n2) is Fourricr transformation of the cor- 
relation function Eq.(3.5) for the momentum difference 
q at the Fermi surface q = pf{^i — ^2) and the sum 
over i, j is implied in Eqs.(3.18), (3.19). The last term, 
Eq.(3.19), corresponding to an effective spin-orbital scat- 
tering describes correlations between the orbital and 
spin scattering. The functional ForbiQn], Eq.(3.17), has 
been previously obtained in Rcf.^". Let us remind that 
Qn = VnAVn in Eqs.(3.17-3.19) is 16 x 16 supermatrix 
with an additional spin structure. 



The functionals, Eqs.(3.17-3.19), have a form typical 
for (T-model in the regime that can be specified as colli- 
sional. From this point of view ni, n2 in Eqs.(3.17-3.19) 
can be considered as momenta of a particle before and 
after collision whereas the integrands determine corre- 
sponding transition probability. Previously we assumed 
that the matrices V'n(r) vary in space sufficiently slowly. 
Now we can clarify this assumption more carefully. 

We notice that its formulation is directly related to the 
conditions under which the coUisional regime described 
by Eqs.(3.17-3.19) should be used. These conditions de- 
pend on the kind of disorder. If disorder is short ranged 
((^-correlated) the scattering may be considered as a colli- 
sion on the length scales exceeding the Fermi wavelength 
Xp. The situation becomes more complicated when the 
disorder is long ranged. This case for the potential dis- 
order was first studied in Refs.^®'^^ and then, with using 
results of these papers and applying the ballistic a-model, 
in Ref.^^. The model with a long ranged magnetic field 
was directly considered in Ref.^^. It was shown that col- 
lisional regime corresponds to lengths exceeding the Lya- 
punov length II = vpX^ , A^ is the rate of divergency of 
the classical trajectories. 

Summing up, the model found in Eqs.(3.17-3.19) 
should be valid for matrices Vn(r) slowly varying over the 
Fermi wavelength Xp in the case of short ranged RMF or 
over the Lyapunov length l^ in the case of a long ranged 
field. 



IV. DIFFUSIVE MODEL. 

In this section we study diffusive limit of the model 
obtained in Eqs.(3.17-3.19) and also calculate the spin 
susceptibility in this limit. 

We start the reduction of the theory given by 
Eqs.(3.17-3.19) to the effective cr-model applicable in the 
diffusive limit noticing that the functionals, Eqs.(3.17- 
3.19), are invariant with respect to rotations U{r) 

g„(r) ^ C/(r)g„(r)C7(r), U{r)U{r) = 1 (4.1) 

provided L'^(r) is the (5-function in the spin space. This 
means that only the charge mode related to fluctuations 
of the electron density remains gapless in the diffusive 
limit. 

The situation changes if one considers a random mag- 
netic field with a fixed direction h, h^ = 1, replacing the 
pair function, Eq.(3.5), by the following one: 

(B(r)B(r'))-2(^)'^r-r'), (4.2) 

where B{r) is the absolute value of the magnetic field. 
Function w{r — r') being considered in 3D should be in- 
dependent of the coordinate along h due to the condi- 
tion diVrB = 0. One can readily see that in this case 
the energies Eqs.(3.17-3.19) are invariant with respect 
to the transformation, Eq.(4.1), if C/(r) commutes with 



the matrix au = (ch) and, consequently, an additional 
soft mode identified with the spin mode along h appears. 
Considering only charge and spin modes, [Qn,crh\ — 0, 
and simplifying Eqs.(3.17-3.f 9) we come to the following 
form of collision term 



Fscatt[Qn\ — — 



iM' 



Sti 



cos 



, 1 + nin2 
1 - nin2 



w(ni - n2)Qm{r)Qn2ir)drdnidn2, 



(4.3) 
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where cos 6 — (ezh), ez is z-ort. Because of the usual 
space symmetry the term Fsp-orh[V], Eq.(3.19), vanishes 
and does not contribute to the energy Eq.(4.3). 

From Eq.(4.3) one can come to the diffusive a model 
singling out the angle modes (5n(r) 

g„(r) = C/(r)g„(r)C7(r), 

that describe the fluctuation of the particle momentum 
n and integrating them out. Calculation of integrals over 
the angle modes Qn is standard and can be found in de- 
tails e.g. in Ref.^". The final expression for the effective 
energy has the usual form 



F[Qir)] = ^Str / [D{V,Q)^ + 2iiuAQ{r)] dv 
where D = VpTtr/2, Ttr is the transport time 



(4.4) 
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)s^6'(l + nin2) + 
xw(ni — n2)(ini(in2 



(1 - nin2) 



(4.5) 



This time determines the relaxation of the angle modes 
and, hence, of the particle momentum n. 

If the magnetic field is (5-correlated in space the trans- 
port time can be explicitly calculated and expressed 
through the time t^^ in the model of spinless electrons 
and perpendicular magnetic field: 



^-1 

'tr 



[cos^ + {g/2f 



(4.6) 



In the case of the magnetic field perpendicular to the 
plane, cos 6 — I, and free electrons, g = 2, the transport 
time Ttr turns out two times smaller than the transport 
time for the spinless electrons. The same result for this 
particular cas has been established in Ref.^^. 

Using Eqs.(3.17-3.I9), (4.4) one can calculate various 
physical quantities. Below we find the spin magnetic sus- 
ceptibility reliing on the diffusive model Eq.(4.4). 

First, we use the Kubo formula and write the mag- 
netic moment due to the spin of electrons as a sum of 
two terms: 



Gf+Jr,r')a.B(r',c.))Gf(r',r)a]dr', 



(4.7) 



Tr corresponds to the trace in the usual spin space. The 
other term niso is given by the same formula taken after 
the replacement (g/4)crB(r', a;) -^ — iA(r',cj)V(.. This 
contribution is related to the both orbital and spin mo- 
tions. It is determined by the angle modes and interac- 
tion between the charge and spin modes. Both factors 
are irrelevant in diffusive limit and the contribution is 
small. 

So, we neglect m^o and calculate nisp. The contribu- 
tion to this quantity given by the second term in Eq.(4.7) 
can be easily found and corresponds to the usual Pauli 
susceptibility. In order to calculate the first term we ex- 
press the product of two Green functions G^ and G'^ 
through the partition function, Eqs.(2.f), (2. II): 



Tr [Gf (ri,r2)^,Gf_Jr2,ri)^.] 



d^Z[a] 
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(4.8) 



where the source a(r, r') is taken as follows: 



a(r,r') = 



0, , «^f ) ) Sir r') 
a2i(r) y ^ ^ 



ai2(r) 



«) 
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Substitution of the partition function Z[a\ taken in diffu- 
sive limit into Eq.(4.8) results in the following expression 
for the susceptibility: 



Xjj(q,w) ^gn\v 



%-2c^^(Tr( (Tf [Q84(q)+g37(q)] 



m(r, w) = msp(r, uj) + mso{v, uj) 
where nisp is a purely spin contribution: 



Tr aJ[Q73(-q)+Q48(-q)] 



(4.10) 



The first term in the brackets is the usual Pauli sus- 
ceptibility. The other term determines the correction to 
this result related to RMF. The angular brackets (■ • •)q 
in Eq.(4.10) imply the averaging of the quantity inside 
them over the matrix Q with the free energy functional, 
Eq.(4.4). 

We calculate this average using the parametrization 



g(r) =K[l + iP{v)\[l-iP{r)\ 



(4.11) 



where the non-diagonal matrix P(r), P(r) = — P(r), 
[P, A] = is the sum of the charge and spin modes 



P(r) = P,h{v) + (Sh)P,p(r) 



(4.12) 



Pc/i(r), Psp(r) are (5- functions in the spin space. Their 
structure should be found in correspondence with 
Eq.(2.14) and condition [Qjfs] ~ 0. This results in the 
following representation: 



Pch{sp) (r) 



B^hisp) (r) 

Bch(sp){r) 



where i3ch(sp)(r) = C'oSj^(sp)(r)Cj, ac/i(sp)(r), ^ch(sp)(r) 
are usual and pi^2ch(sp) ~ Grassman matrices 4x4 unit, 
in spin and particle-hole spaces: 



ach(sp)i'r) 



ach(sp)iy^) 






^ch{sp] 



(r) 



bcHsp){r) - [ 6:,(,^)(r) 



/ N / Pl,2ch(sp)(r) \ fA TA^ 
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In Gauss approximation the charge and spin modes 
give independent contributions to the energy Eq.(4.4). 
Standard calculation of the average in Eq.(4.10) in the 
same approximation gives: 

Xij(q,w) ^ gp%v{5ij - hih-j) + gp\v-—^-——hihj 

(4.15) 

In the particular case when RMF has only the perpen- 
dicular component, this expression is in agreement with 
the result of Ref^^ . 



V. DISCUSSION 



In the present work we considered a two-dimensional 
electron gas placed in a nonuniform (random) magnetic 
field (RMF) using method of superbosonization Rcf.^^. 
The problem was studied in a quite general formulation 
with taking into account interaction of the magnetic field 
with the spin of electrons. The direction of the magnetic 
field was assumed random in space and the value of g- 
factor - arbitrary. 

In the present paper we used a new method proposed 
in Rcf.^^. This method is based on the exact transforma- 
tion of the initial field theory written in terms of a func- 
tional integral over supervectors -0 to a theory described 
by a functional integral over supermatrices Q and, by the 
analogy with quantum field theory, can be called super- 
bosonization. As the method is exact it enabled us to 
deal with both the short and long range disorder. 

The found theory turned out to be invariant with re- 
spect to the rotations in the superspace Q -^ VQV, V 
are unitary supermatrices VV — 1, provided that they 
commute with Hamiltonian: [-ff , y] = 0. This means 
that the low energy limit of the theory is described by 
the rotations V with a small commutator [//, y]. Con- 
ditions which should be imposed on the supermatrices 
V to provide the smallncss of the commutator depended 
on the character of the terms contained in Hamiltonian 
H , e.g. on the kind of disorder. Singling out the massive 
modes related to the fluctuations of the advanced and re- 
tarded blocks q in the representation Q ~ VqV where q - 
arbitrary diagonal supermatrices [q, A] = and integrat- 
ing them out we came to an effective energy functional 
in terms of the rotations V only and obtained some non- 
linear a model. 

The derivation of the a- model was carried out for an 
arbitrary Hamiltonian H and did not depend on its ex- 
plicit form. Particularly, it might be used for both short 
and long range disorder. Explicit form of the model as 
well as the limits of its applicability depended on charac- 
ter of terms contained in Hamiltonian H . In this paper 
we dealt only with a weak RMF without an additional po- 
tential disorder and considered the theory in the regime 
that may be called coUisional. This regime was achieved 
differently depending on correlation length of RMF. If the 
field was (5-correlated the model was valid at distances 
larger than wavelength Xp- In the case of long range 
RMF they should be larger than the Lyapunov length 
II- 

After averaging the free energy functional over RMF 
we obtained the cr- model with a collision term consist- 
ing of three contributions that could be called orbital, 
spin and spin-orbital. The first contribution came from 
the part of the Hamiltonian acting on the orbital motion, 
the second one- from the part acting on spin, whereas the 
last one came in the result of joint averaging of them and 
describes correlations between the spin and orbit scatter- 
ing. 



The coUisional term turned out to be invariant with 
respect to the rotation Qn -^ UQU provided that U was 
i5-function in the spin space. This allowed us to conclude 
that in general only the mode related to the charge trans- 
port or charge mode should survive in the diffusive limit. 
At the same time, we have shown that in the particu- 
lar case of RMF with a fixed direction h one more spin 
mode along should be regarded as well. We obtained the 
(T-model in the standard diffusive form and found the 
transport time Ttr- Due to the usual space symmetry the 
spin-orbital collision term vanished in the model with a 
not fluctuating fixed RMF and did not contribute to Ttr- 

Finally, in diffusive limit, we calculated the spin 
magnetic susceptibility. The part of the susceptibility 
transversal to the direction h is given by the usual Pauli 
expression whereas the longitudinal one has an additional 
contribution proportional to the diffusion propagator. 

We acknowledge a support from SFB 491 and 
SFB/Transregio 12. 
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